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Consequences of anisotropy (variation of orbital occupation) and magnetism, and their coupling, 
are analyzed for LSDA+U functionals, both the commonly used ones as well as less commonly 
applied functionals. After reviewing and extending some earlier observations for an isotropic inter- 
action, the anisotropics are examined more fully and related to use with the local density (LDA) 
or local spin density (LSDA) approximations. The total energies of all possible integer configura- 
tions of an open / shell are presented for three functionals, where some differences are found to be 
dramatic. Differences in how the commonly used "around mean field" (AMF) and "fully localized 
limit" (FLL) functionals perform are traced to such differences. The LSDA+U interaction term, 
applied self-consistently, usually enhances spin magnetic moments and orbital polarization, and the 
double-counting terms of both functionals provide an opposing, moderating tendency ( "suppressing 
the magnetic moment"). The AMF double counting term gives magnetic states a significantly larger 
energy penalty than does the FLL counterpart. 

PACS numbers: 



I. INTRODUCTION 



Density functional theory (DFT) and its associ- 
ated local (spin) density approximation [L(S)DA] is 
used widely to describe the properties of a wide va- 
riety of materials, often with great success. However 
there exists a class of materials which are poorly de- 
scribed, sometimes qualitatively, by LDA. These so- 
called strongly correlated materials typically contain 
atoms with open d or / shells, in which the corre- 
sponding orbitals are in some sense localized. The 
LSDA+U approach was introduced by Anisimov, 
Zaanen, and Andersen^ to treat correlated materials 
as a modification of LDA ('on top of LDA') that adds 
an intra-atomic Hubbard U repulsion term in the 
energy functional. Treated in a self-consistent mean 
field ('Hartree-Fock') manner, in quite a large num- 
ber of cases the LDA+U result provides a greatly 
improved description of strongly correlated materi- 
als. 

At the most basic level, the LSDA+U correction 
tends to drive the correlated orbital m occupation 
numbers n m<7 (er denotes spin projection) to integer 
values or 1. This in turn produces, under appro- 
priate conditions, insulating states out of conduct- 
ing LSDA states, and the Mott insulating state of 
several systems is regarded as being well described 
by LSDA+U at the band theory level. Dudarev et 
al£ and Petukhov et a/4 provided some description 
of the effect of the spin dependence of two different 
double counting terms within an isotropic approx- 
imation. Beyond this important but simple effect, 
there is freedom in which of the spin-orbitals (ma) 



will be occupied, which can affect the result consid- 
erably and therefore makes it important to under- 
stand the effects of anisotropy and spin polarization 
in LSDA+U. After the successes of providing realis- 
tic pictures of the Mott insulating state in La2CuC>4 
and the transition metal monoxidesji the anisotropy 
contained in the LSDA+U method produced the 
correct orbitally ordered magnetic arrangement for 
KCUF3 that provided an understanding of its mag- 
netic behavior A 

The anisotropy of the interaction, and its con- 
nection to the level of spin polarization, is a topic 
that is gaining interest and importance. One exam- 
ple is in the LSDA+U description of the zero tem- 
perature Mott transition under pressure in the clas- 
sic Mott insulator MnO. The first transition under 
pressure is predicted to be£ an insulator-insulator 
(not insulator-metal) transition, with a S = I — > 
S = h moment collapse and a volume collapse. The 
insulator-to-insulator aspect is surprising, but more 
surprising is the form of moment collapse: each or- 
bital remains singly occupied beyond the transition, 
but the spins of electrons in two of the orbitals 
have flipped direction. This type of moment col- 
lapse is totally unanticipated (and hence disbelieved 
by some), but it is robust against crystal structure 
(occurring in both rocksalt and NiAs structures) 
and against reasonable variation of the interaction 
strength. Detailed analysis indicates it is a product 
of the anisotropy of the LSDA+U interaction and 
the symmetry lowering due to antiferromagnetic or- 
der. 

Another unanticipated result was obtained^ in 



LaNiC>2, which is a metal experimentally. This com- 
pound is also a metal in LSDA+U over a very large 
range of interaction strength U, rather than revert- 
ing to a Mott insulating Ni 1+ system which would 
be isovalent with CaCuC>2. For values of U in the 
range expected to be appropriate for the Ni ion in 
this oxide, the magnetic system consists of an atomic 
singlet consisting of antialigned d x i_ y i and d z 2 spins 
on each Ni ion. Again the anisotropy of the interac- 
tion evidently plays a crucial role in the result, with 
its effect being coupled thoroughly with band mixing 
effects. 

The addition of a Hubbard U interaction also in- 
troduces the need for "double counting" correction 
terms in the energy functional, to account for the 
fact that the Coulomb energy is already included (al- 
beit more approximately) in the LSDA functional. 
All double counting schemes subtract an averaged 
energy for the occupation of a selected reference 
state depending only on {N a }, which largely can- 
cels the isotropic interaction of the Ei term Eq. @ . 
Several forms for these double-counting terms have 
been proposed jiiLS but primarily two are commonly 
used. These LDA+U functionals are most often re- 
ferred to as around mean field (AMF) and the fully 
localized limit (FLL), which is also referred to as the 
atomic limit (AL). The distinctions between these 
forms have attracted some discussion, but without 
consideration of the full anisotropy of the interac- 
tion. 

The need for double-counting corrections is not 
unique to the LDA+U method; any other method 
that adds correlation terms to the LSDA func- 
tional, such as the dynamical LDA+DMFT (dynam- 
ical mean field theory) approach, also will require 
double-counting corrections. This is an unfortunate 
consequence of LDA's success; LDA works too well, 
even in correlated systems where it usually gets in- 
teratomic charge balance reasonably, to just throw 
it away4£ The common approach has been to use 
LSDA for correlated materials and include a double- 
counting correction. There are techniques being de- 
veloped which do not build on a correction to DFT- 
LDA, but it remains to be seen whether these ap- 
proaches will be successfully applied to a broad range 
of solid-state materials. 

Although there has been much study on the per- 
formance of these LDA+U functionals in the context 
of real materials, and an early review of the method 
and some applications was provided by Anisimov, 
Aryasetiawan, and Lichtenstein^ relatively little has 
been done to understand, qualitatively and semi- 
quantitatively, how the functionals operate based 
solely on their energetics, distinct from DFT-LSDA 
effects. In this paper we analyze the functionals that 
are commonly used, as well as others which were in- 
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troduced early on but are not so commonly used. 
Some of the nomenclature in the literature is con- 
fusing, so we try to clarify these confusions where 
we can. 

II. THE LSDA+U CORRECTION AE 

The LDA+U functional is usually coded in a form 
in which the choice of coordinate system is irrele- 
vant, often referred to as the rotationally- invariant 
form. 4 This form involves Coulomb matrix elements 
that have four orbital indices, and the orbital occu- 
pation numbers are matrices in orbital space (viz. 
n mm ')- One can always (after the fact) rotate into 
the orbital Hilbcrt space in which the occupations 
are diagonal, in which case the interactions have only 
two indices. In our discussion we will work in the di- 
agonal representation. 

The LDA+U functionals considered here can all 
be written in the form 

AE = Ej - E dc , (1) 

where the direct interaction is 

Ei = \ W ^Cn m ,n m , a , (2) 

and Edc is the double-counting correction. The 
Coulomb matrix elements are given in terms of the 
direct and (spin-dependent) exchange contributions 

as 

^mm' = {Umm f Jmm' $ a^a'*) • (3) 

By the convention chosen here, Ei and Ed c are both 
positive quantities as long as the constants U and J 
(which define the matrix elements U mm > and J m m' 
but are not the same) as chosen conventionally, with 
U much larger than J. 

Note that the orbital+spin diagonal term has been 
omitted in Eq. [2] - there is no self-interaction in 
Ej. However, it is formally allowed to include the 
diagonal 'self-interaction' term, because the matrix 
element vanishes identically (self-interaction equals 
self-exchange: U mm = J mm ), and it can simplify 
expressions (sometimes at a cost in clarity) if this 
is done. The double counting correction depends 
only on the orbital sum N a , which appears up to 
quadratic order. A consequence is that it will con- 
tain terms in n ma n mtTl which are self-interactions. 
Thus while the LSDA+U method was not intended 
as a self-interaction correction method, it is not 
totally self-interaction free. In fact, the underly- 
ing LSDA method also contains self-interaction, and 
the double-counting term may serve to compensate 
somewhat this unwanted effect. We discuss self- 
interaction at selected points in this paper. 
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LDA+U 
Functional 


E dc = E dc DFTXC 

(rewritten) Functional 


Fl-nS 


frtf-mfttf = - HffiL^fi)* LDA 


Fl-S (AMF) 


LUN^^LLL^m = LUN^^L^m LSDA 


FLL 


±UN(N- 1) - hJY. a N a {N a - 1) = ±UN(N -l)-%JY,ANt-N a ) LSDA 


FLL-nS 


\UN{N - 1) - \JN(N - 2) = i[/iV(iV - 1) - |JE CT ((f ) 2 - W CT ) LDA 



TABLE I: The double-counting terms of various LDA+U functionals. In the second expression two of them are 
rewritten to reflect how they are (somewhat deceptively) identical in form, but in one case a distinction between 
spin-up and spin-down (relative to half of N: N a <-> N/2) is made. Note that while the first two forms appear to 
contain an isotropic self-interaction [^UN 2 rather than ^N(N — 1)] they are derived from a form which explicitly has 
no self-interaction between the orbital fluctuations 5 Ma- See text for more discussion. 



A. Short formal background to the LSDA+U 
method. 

The "LSDA+U method" is actually a class of 
functionals. Each functional has the same form of 
interaction Ej, with differences specified by 

(1) choice of the form of double counting term. 

(2) choice of constants U and J. For a given func- 
tional, these are 'universal' constants like ft, m, e, i.e. 
they are not functional of the density in current im- 
plementations. Possibilities for doing so, that is, de- 
termining them self-consistently within the theory, 
have been proposed^ 

(3) choice of projection method to determine the 
occupation matrices from the Kohn-Sham orbitals. 
Given identical choices for (1) and (2) above, there 
will be some (typically small) differences in results 
from different codes due to the projection method. 

The occupation numbers (or, more generally, ma- 
trices) are functionals of the density, n mrT [p] , through 
their dependence on the Kohn-Sham orbitals. Then, 
whereas in LSDA one uses the functional derivative 

LSDA: dE ^DA[{Ps}] (4) 

in minimizing the functional, in LSDA+U the ex- 
pression generalizes to 

LSDA + U: 9(E LSDA{Ps} ] + AE [{ n ms[Ps]}] ) 

Since the resulting spin densities p s are changed by 
including the AE correction, the change in energy 
involves not only AE but also the change in Els da- 
In practice, there is no reason to compare E LS da+u 
with Elsda as they are such different functionals. 
However, in the following we will be assessing the im- 
portance of the choice of the double counting term 
in the LSDA+U functional, and it is of interest to 
compare, for fixed U and J, the energy differences 
between LSDA+U functionals differing only in their 
double counting terms in order to understand the 



differing results. Even if the set of occupation num- 
bers turn out to be the same (a situation we consider 
below) , the densities p a will be different and the dif- 
ferences in Elsda may become important. 

As with the non- kinetic energy terms in Elsda, 
the functional derivatives of AE lead to potentials 
in the Kohn-Sham equation. These are non-local 
potentials, which (via the same projection used to 
define the occupation numbers) give rise to orbital- 
dependent (nonlocal) potentials 

._ dAE _ j _ dc 

Vmo- — rv — V 7na v mcn 

v La = Cm'»mV. (6) 

ra' a' ^raa 

The double counting orbital potential is discussed 
later. 

The corresponding contribution to the eigenvalue 
sum E sum is 

AE su7n ^ ^mff^mdi (7) 

which is subtracted from the eigenvalue sum to ob- 
tain the Kohn-Sham kinetic energy. However, there 
are indirect effects of the orbital potentials that af- 
fect all of the kinetic and (LSDA) potential energies; 
these will be different for different AE functionals 
because the orbital potentials, which depend on the 
derivative of AE and not simply on the values of 
n ma , differ for each functional. This makes it neces- 
sary, for understanding the effects of the AE correc- 
tion and the change in energy, to analyze the orbital 
potentials. We provide a brief discussion in Sec. V. 

B. Fluctuation forms of LSDA+U 

First we consider the class of functionals that can 
be written in what is termed here as a fluctuation 
form. The original LDA+U functional was intro- 
duced in 1991 by Anisimov, Zaanen and Andersen^ 



and was written as 

AE Fl ~ nS = 1 £ W°C(nm*-n)(n m , a ,-n), 

(8) 

where n — N co " /2(2l + 1) is the average occupa- 
tion of the correlated orbitals. (Henceforth N = 
N COT1 '.) Note that the energy is changed only ac- 
cording to angular 'fluctuation' away from the (spin- 
independent) angular average occupation. This form 
is properly used with LDA (the 'LDA averages' n are 
the reference) and not LSDA. This form was origi- 
nally advocated with generic (U — JS^^i) matrix el- 
ements instead of the full Coulomb matrix, but we 
use the full W™ m , here for comparison with other 
functionals. 

In 1994, Czyzyk and Sawatzky^ introduced a 
change to |5]) and also proposed a new functional. 
The motivation for changing {SJ was to use an LSDA 
exchange-correlation functional to treat spin split- 
ting effects rather than LDA. This change motivated 
the following equation, 

AE Fl - s = - J2 W™ m (n ma -n a )(n m , a , - n a .) 

ma^ra' a' 

= AE AMF (9) 

where fi a — N a /(2l + 1) is the average occupation of 
a single spin of the correlated orbitals. Here the en- 
ergy correction is due to angular fluctuations away 
from the spin-dependent angular mean, and hence 
must be used with LSDA. We point out that the au- 
thors in Ref. [| refer to Eq. © as ^lda+amf 
and Eq. © as #lsda+amf^ Thig wording may 

have caused subsequent confusion, due to the way 
these terms have come to be used, and also be- 
cause a discussion of the "+U" functionals requires 
explicit specification of whether LDA or LSDA is 
being used just to understand which functional is 
being discussed. Also confusing is that Solovyev, 
Dederichs, and Anisimov-ii rejustified Eq. [5] using 
"atomic limit" terminology. 

The fluctuation forms of LSDA+U are automati- 
cally particle-hole symmetric, since n ma — > 1 — n m<T , 
rifj y 1 tl gives n ma — n a ^ —(n rn<T — h a ) and the 
expression is quadratic in these fluctuations. The 
general form of Eq. [T]need not be particle-hole sym- 
metric. 

Many authors (present authors included) have 
used the term LDA+U where the term LSDA+U 
would be more appropriate, which is especially con- 
fusing when discussing the AMF functional. We 
choose to depart from this confusing nomenclature 
by giving (jHJl and (© unique names specifying their 
fluctuation forms, and their connection to LSDA (Fl- 
S) or to LDA (Fl-nS). We collect the double-counting 
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terms for the various functionals, along with their 
connection to LDA or LSDA, in Table I. 



C. The Fully Localized Limit (FLL) Functional 

The second functional introduced by Czyzyk and 
Sawatzky® is the FLL functional. (A J=0 version 
of FLL was introduced in 1993 by Anisimov et alX) 
The authors referred to it (confusingly, as terminol- 
ogy has progressed) as the "around mean field" func- 
tional but the atomic limit double counting term; in 
the literature it is now referred to as the atomic limit, 
or fully localized limit (FLL) functional. This func- 
tional cannot be written in the fluctuation form of 
the previous two functionals (the fluctuation form 
is exhausted by the -S and -nS cases). The FLL 
functional is written in the form of (|TJ) , with the 
double-counting term given in Table |TJ 

There is yet another LDA+U functional that is 
available, which was introduced in 1993 by Anisi- 
mov et alX There is no clear name for it, but since 
it can be obtained by using N a — N/2 in Ed c for 
FLL, one might consistently refer to it as FLL-nS, 
corresponding to FLL with no spin dependence. The 
authors in Ref. 7 indicate that this functional is to 
be used with LDA, in accordance with the lack of 
spin dependence in the double counting term. 



D. Implementation of LSDA+U 
in Some Widely Used Codes 

The Fl-nS functional is implemented in the 
Wien2k code, as nldau = 2, and called HMF 
(Hubbard in Mean Field) ^ however it is appar- 
ently not often used. 

The Fl-S (AMF) functional is implemented in the 
Wien2k codei 2 . as nldau = and the FPLO codei^ 
as AMF. It is also available in the Abinit code 14 
when using a PAW basis seti^ by setting usepawu 
= 2. 

The FLL functional is implemented in several 
general-purpose DFT codes, such as Wien2k 
(nldau = l) 12 , FPLO (select AL in fedit)^, VASP, 
PW/SCF, and Abinit (usepawu = 1) when using a 
PAW basis set. 

The FLL-nS functional is available in VASP. 



E. General Remarks 

When the Fl-nS and Fl-S functionals are writ- 
ten in their fluctuation form, there is no separate 
double-counting term, hence one does not need the 



double-counting interpretation. They can of course 
be expanded to be written in the 'interaction mi- 
nus double-counting' form of (J2]), which is useful 
especially for comparing with functionals that can 
only be written in that form. A comparison of the 
double-counting terms is given in Table |TJ Reducing 
all to interaction minus double-counting form makes 
the difference between the functionals most evident; 
since they all have the same "direct-interaction" 
term, the only difference between the functionals is 
what double-counting energy is used; the uninterest- 
ing tail seems to be wagging the exciting dog, which 
is in fact the case. The double-counting terms can 
be reduced to dependence only on N and N a thanks 
to summation rules (there is at least one free index) 
on the U mm i and J mm ' matrices, 



With J = 0, the fluctuation functionals simplify to 



m 

2_j Jn 



(2l + l)U 
U + (21) J, 



(10) 
(11) 



that is, the sum over any column (or row) of the U 
and J matrices is a fixed simple value, which depends 
on the input parameters U and J. One can then 
simply see that a sum over a column of W is (2l+l)U 
if a j= a' and 21 (U - J) if a = a'. 

The U mm i and J mm i matrices satisfy, by defini- 
tion, U mm — J mm = 0, so that there is no self- 
interaction, whether or not the (vanishing) diagonal 
term ma = m'a' is included in the interaction term. 
As mentioned earlier, the following analysis assumes 
the the occupation matrix has been diagonalized. 
While this can always be done, the transformed ma- 
trix elements U m m' and J m m' will not be exactly 
what we have used in Sec. IVH 



III. ANALYSIS OF THE FUNCTIONALS 



1. The J — simplification. 



AE 



Fl-h 
J=0 



u 

~2 

U_ 
~2 



E 5fl r> 

^ Sn mc 



42> 



- y^(fa mCT ) 

rn a 

= -^r 2 <o, (12) 



because the sum of fluctuations vanishes by defini- 
tion for either form n — nS or S; note the 'sign 
change' of this expression when the diagonal terms 
are added, and subtracted, to simplify the expres- 
sion. Here T 2 is the sum of the squares of the fluc- 
tuations, bounded by < r 2 < N. For integer oc- 
cupations the energy corrections for Fl-nS and Fl-S 
(AMF) can be written 



AE 



Fl-S 
J=0 



U 
~2 

U 



N(l 



M 2 



2(21 + 1) 



A25j:=o = --N(l-n). 
There are two things to note here. 



(13) 



1. In Fl-nS, the energy is independent of both 
the spin and orbital polarization of the state, 
which lacks the basic objective of what 
LSDA+U is intended to model. Considering 
the form of its double counting term (see Table 
|T| with its self-interaction term (proportional- 
ity to ^V 2 ), Fl-nS for J=0 becomes simply a 
self-interaction correction method. 

2. In Fl-S (AMF), configurations with magnetic 
moments are energetically penalized, propor- 
tionally to U and quadratically with M. In 
later sections we will discuss the partial can- 
cellation with the LSDA magnetic energy. 

Under the same conditions, the FLL functional 
becomes 

AE FLL = ^Yj ""-t 1 ~ ^ °- ( 14 ) 



It is not uncommon for practitioners to use 'effec- 
tive' values U = U — J, J = and insert these con- 
stants (for U, J) into LDSA+U. For J = 0, of course 
Hund's coupling (intra-atomic exchange) is lost, but 
J also controls the anisotropy of the interaction, and 
for J = anisotropy also is lost {Umm 1 = U as well 
as Jmm' = for m ^ m'). This case is relatively 
simple, it seems it should provide the "big picture" 
of what LSDA+U does with simple Coulomb repul- 
sion, and it has been discussed several times before. 



Solovyev et al^ noted the important and easily rec- 
ognizable characteristics of this expression. Besides 
being non-negative, for integer occupations the en- 
ergy vanishes. It is a simple inverted parabola as 
a function of each n mCT . From the derivative, the 
orbital potentials are linear functions of n mCT , with 
a discontinuity of U when n ma crosses an integer 
value. These characteristics underlie the most basic 
properties of the LSDA+U method: integer occupa- 
tions are energetically preferred, and discontinuities 



in the potentials model realistically the Mott insula- 
tor gap that occurs in strongly interacting systems 
at (and only at) integer filling. 



2. J 7^ 0, but Isotropic 

Simplification of the full expression for a func- 
tional results by separating out the isotropic parts 
of the interaction: 



Un 



U + AU„ 



J rr 



U6 

III in + J(l - 6 mm >) + AJ n 



(15a) 



(15b) 



The isotropic parts simplify, giving 

mer 



U — J 



Nn + AE a 



(16) 



AE 



Fl-S 



4 21 + 1 



U -J 



Nn + AE a 



(17) 



ae fll = _R_±y n 2 a + E—L N 

2 / j Trier 2 
mer 



+AE aniso (18) 
with the universal anisotropy contribution 

As aniso = 1 AW'Cn^n^,, (19) 

mm'erer' 

AW™ m , = AU mm , -AJ mm ,8„„,. (20) 

is the anisotropic part of the interaction matrix ele- 
ments. These equations, up to the AW term, are the 
extensions of Eq. (fl~3| to include isotropic exchange 
in explicit form. 

The first term in each of these expressions contains 
— ^Un^iU = U — J) and hence has the appearance 
of a self-interaction correction. Since the diagonal 
term of the interaction Ei is specifically excluded, 
it does not actually contain any self-interaction; in 
fact, the sign of the interaction Ej is positive. (The 
double counting term does contain terms quadratic 
in N which must be interpreted as self-interaction.) 
Nevertheless, the rewriting of the functional leads 
to a self-interaction-like form, and that part of the 
functional will have an effect related to what appears 
in the self-interaction-corrected LDA method, but 
by an amount proportional to U rather than a direct 
Coulomb integral, and depending on the difference 
of n ma from the reference occupation, see Sec. IV. 



3. Fl-nS 

For Fl-nS, if we are restricted to integer occupa- 
tions (so n 2 na = n ma ), then T 2 depends only on N, 
so the first term in AE Fl ~ nS above depends only on 
N. Then, up to corrections in AU and A J, the state 
with the largest total spin moment will be favored; 
this is Hund's first rule. In fact, even with the AU 
and AJ terms, the — JM 2 /4 term is still strongly 
dominant. Except for N = 7, there are many ways 
to arrange electrons in orbitals which maximizes S. 
Energy differences between these arrangements arise 
only from anisotropy (AU and AJ) and spin-orbit 
coupling. 



4. Fl-S 

In Fl-S, instead of having the — JM 2 /4 term 
from Fl-nS which favors magnetism, there is a term 
M2i+i) w hich opposes magnetism. This term 
(as in the J = case) comes from the occupation 
variance which wants to evenly distribute electrons 
across both spin channels. Within LSDA there is 
something like a Stoner term of the form — jIM 2 
which will compete with this Fl-S magnetic penalty. 
We return to this aspect in later sections and the 
appendix. 



5. Spin-orbit Coupling; Particle-Hole Symmetry 

Without spin-orbit interaction, for a given N 
there are many states that are degenerate for both 
double counting schemes. Every value of TV has at 
least four degeneracies, those with ±L Z ,±S Z . 

Any state which has the same number of spin up as 
spin down electrons [M = 0) gives the same energy 
from Fl-nS and Fl-S, since then nj = = n (the or- 
bital potentials are distinct, however). Of course this 
fixed N, M=0 specification may contain many dif- 
ferent configurations. Looking at results mentioned 
later, for Fl-S the ground state for an even number 
of electrons is S z = (so n a = n), thus the con- 
figuration which gives the Fl-S ground state has the 
same energy in Fl-S and Fl-nS. 



IV. FRACTIONAL OCCUPATIONS 

Here we briefly discuss the effect of non-integer oc- 
cupations in LSDA+U. Taking a general set of occu- 
pations as {n m(T }, we define a set of integer occupa- 
tions, {n ma }, and the fractional part of the occupa- 
tions as 7 TO(T = n ma —h ma . For illustration purposes 



we will choose the simplest possible scenario, where 
charge is transferee! to an empty orbital a from an 
occupied orbital 6 both of the same spin, so that 
< 7a| — — 7&T' "at ~ and = 1. With this 
selection, N a is unchanged (and therefore, N and M 
as well) so that Ed c is unchanged. Thus, the effect 
of the charge transfer is entirely contained in the Ej 
term. Expanding E\ for the general occupation set 
gives 

Ei[{n m<J }] - £ 7 [{n mCT }] = U la] {l - 7aT ) (21) 
for the J — case, and for J ^ we find 
Ei[{n ma }} - £/[{n mcr }] 



ab la] 



(22) 



The dominant term in 



is where ma = b f. 

■V, 



This term gives a contribution W'. 7 a f ~ U^ a ^ 
(since U >> J for typical parameter choices, 
where other terms give contributions proportional to 
- W b T ^) 7aT oc J 7aT . The term with ma = a t 
is killed off by the factor of n a f, and the term in j 2 
is significantly smaller than the others for 7 < 0.5. 

This shows that there is an energy penalty for 
fractional occupation, proportional to U and lin- 
ear in 7 at small 7. Thus, in configuration space, 
the LSDA+U functionals have many local min- 
ima around configurations with integer occupations. 
This result is fairly general. Even for charge trans- 
fer between orbitals of opposite spins, the linear 
energy penalty in 7 is still dominant over any ad- 
ditional terms coming from the double-counting or 
spin-orbit. 

In practice, this gives the possibility that 
LSDA+U will get 'stuck' in a local minimum with 
some configuration that may not be the true ground 
state. This behavior is not uncommon; LSDA+U 
has been reported^ to find multiple local minima 
depending on the starting configuration. 



V. ORBITAL POTENTIAL MATRIX 
ELEMENTS 

Up to now only the energy functionals themselves 
were discussed. Now we return to the derivatives, 
the orbital potentials v ma . It is simple to derive 
the exact expressions, and the interaction term Ej 
common to all forms gives a potential Av ma - which 
depends only on the occupations of the other orbitals 
n m i a i, m'a' ^ ma. The potential resulting from the 
double counting term is functional specific, and may 
contain a contribution from n ma itself, i.e. a self- 
interaction. 
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We confine our observations here to the subdivi- 
sion (introduced just above) of the interaction into 
an unitarily invariant isotropic part, and into an 
anisotropic part Eq. (|15|) that is much smaller and 
more difficult to analyze. As for the energy itself, it 
is convenient to add and subtract the diagonal self- 
Coulomb and self-exchange, which makes the effect 
of the potential much more transparent at the cost 
of introducing the misleading self-interaction inter- 
pretation. 

The potential matrix elements are 



Av 



Fl-nS 
ma 



Av, 



Fl-S 



-(U - J) [n ma - n] - -Ma 
+Av%r, (23) 

,TT T\ F - 1 U-JM 

-{U -J) [n ma - n a \ + 2 2T+1 CT 

(24) 



am so 
ma 1 



Av. 



FLL 



+Av 
-(U-J) 



Avl 



with the anisotropic potential term 



Av 



aniso 
ma 



E 



(25) 



(26) 



The main occupation number dependent term, pro- 
portional to n ma , has a self-interaction appearance 
and effect, as discussed above for the functionals. 
The differences in this term arise from the "refer- 
ence" occupation with which n m<7 is compared to 
determine the potential shift. The "fluctuation" 
rima — Tiret is smaller for Fl — S (AMF) than for 
Fl — nS because the occupation for a given spin di- 
rection tends to be closer to n a than to n. The 
reference occupation for FLL is, like Fl — nS, spin- 
independent, in fact, the reference is half-filling. In 
this sense, FLL seems more like a single-band Hub- 
bard model treatment than the other two function- 
als. 

The other difference that is evident in this form is 
the spin dependence. Fl — nS additionally has a spin 
orientation dependent potential shift proportional to 
J and to M (similar to an LSDA treatment, but us- 
ing J instead of the Stoner /) and enhances spin- 
splitting of the eigenenergies e accordingly. In Fl- 
S (AMF) the analogous term is + (U - J) 2 (2i+i) u i 
with a sign that impedes magnetism. It can be sim- 
plified to w f Mct when U « 2 (I + 1) J . This ex- 
pression illuminates the reason that AFM is some- 
times found to decrease the magnetic moment: this 
term more or less cancels the spin splitting of LSDA 
due to the opposite sign. What is left is a splitting of 
occupied and unoccupied levels due to the n ma term, 
which is almost independent of M. The effect is to 
support a spin-polarized solution, but provide little 




N 



FIG. 1: (Color online) Angular momentum values of 
S Z ,L Z , and J z of the lowest energy state for (a) AMF 
(Fl-S) and (b) FLL, with spin-orbit coupling. Parame- 
ter values are U = 8, J = 1, I = 0.75. The AMF (Fl-S) 
curves do not follow Hund's rules, because the Stoner pa- 
rameter is too small. FLL follows Hund's rules exactly 
with these parameters. 
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FIG. 2: (Color online) Shown here is AE FLL plotted vs. 
AE Fl ~ s for each of the 3432 configurations of N = 7 
electrons, using U — 8, J — 1,1 — 0.75, all in eV. The 
ordering of states is shown for Fl-S by counting from left 
to right, and for FLL by counting from bottom to top. 
Open squares show values for U — 7 and J = 0. 



discrimination between different M. Since the spin 
polarization energy does not favor large M, we end 
up with a tendency of a near degeneracy of different 
M values, as we already pointed out from purely en- 
ergetic arguments. For the case of a half-filled fully 
polarized shell n ms = 8 a ,i (the case N = M = 7 in 
Section |VT|) the potential matrix vanishes, which can 
be seen from n — -| 2 (2i+i) = ^ = 5" However, at 
the same time the energy contribution also vanishes 
AE Fl ~ s = (for integer occupations) and the Fl-S 
functional has no effect at all. 

The SIC term in FLL splits occupied and unoc- 
cupied states symmetrically, while in the fluctuation 
functionals the splitting happens with respect to the 
averaged occupation, which is seen in the overall en- 
ergy positions in Fig [3J 



VI. NUMERICAL RESULTS 

Following common terminology, for the remainder 
of the paper we refer to the Fl-S functional simply 
as the AMF form. We have taken values for U mm i 
and Jmm' (used for Eu) from Ref. IT?] (recalculated, 
to include more significant figures). These matrices 
are generated using U = 8 and J — 1 (values typical 
of rare earths) following the procedure given in the 
appendix of Ref. Il| 

In our analysis of the AMF and FLL functionals, 
which are based on an LSDA reference state, we in- 
clude a Stoner term 

E{M) = -jIM 2 (27) 



to model the magnetic effects of LSDA on the total 
energy. The addition of this term helps to give a 
picture of the degree to which the functionals repro- 
duce Hund's first rule. Typical values oil for ionized 
lanthanidcs are 0.75 eV, so we use this value for the 
calculations of this section. Further discussion of the 
Stoner I is included in the appendix. 

Spin-orbit interaction is included in the form 

E so = \S-L^J2 S * L *> ( 28 ) 

I I I XT 

where the second form applies when only z- 
components of moments are treated, as is done in 
current implementations of the LSDA+U method. 
Due to this restriction, LSDA+U often does not pro- 
duce the correct multiplet energies in the atomic 
limit. The visible result in LSDA+U band struc- 
tures is splittings of occupied, or unoccupied, corre- 
lated suborbitals that can be as large as a few times 
J, and understanding the splittings is not straight- 
forward. For 4/ systems these splitting a 17 ! 18 may 
not be of much interest unless one of the correlated 
bands approaches the Fermi level. In heavy fermion 
compounds, for example, LSDA+U results are used 
to infer which parts of the Fermi surface has a larger 
amount of / character^ The same effects (eigen- 
value splittings) occur in 3d or 5/ systems, however, 
where they are expected to become more relevant 
but are masked by stronger banding tendencies. 

Here we consider values of A of and 0.2 eV. The 
magnitude of the spin-orbit interaction is not critical 
to the results; it mainly serves to break degeneracies. 
Without the spin-orbit interaction, the ground state 



for any of the functionals at a given N is degenerate 
with several other states. For instance with N = 6, 
the AMF functional has states with L z = 1, S z = 
and L = 11, S z = with the same lowest energy. 

In Fig. [1] are the ground states for both AMF 
and FLL with U = 8, J = 1 and I = 0.75. The 
FLL and Fl-nS (not shown) schemes both reproduce 
Hund's rules exactly with these parameters. AMF 
does not reproduce Hund's rules (in fact penalizes 
magnetism) until / is increased to around 1.5, which 
is somewhat larger than reasonable values of I. If 
one expects LSDA+U to reproduce Hund's rules, 
then the AMF scheme performs rather poorly. For 
instance, at N = 7, Hund's rules ask that all elec- 
trons be spin-aligned, but the AMF ground state 
has only one unpaired spin due to the magnetic 
penalty appearing in Eq. (|13[) . With these param- 
eter choices, U/(2l + 1) > J, so the AMF magnetic 
penalty wins over the Stoner energy. This is likely 
to be the case for 3d transition metals as well, since 
U-sdj (2Z + 1) ~ leV, but it may not be as significant 
since I for 3d elements is larger. 

We examine the energetics in more detail in Fig. 
H where AE for the AMF and FLL functional is 
plotted for every configuration for N — 7. The con- 
figurations fall into separate lines for each spin mo- 
ment M, since E& c depends only on iV and M for 
both functionals. For the case of J = 0, all the 
states with a particular M value collapse to a single 
energy value (the orbital index loses any impact), 
this is shown with the open squares. A value of I 
was chosen so that the cancellation discussed in the 
previous paragraph is slightly broken. 

If we examine the J = case first (the large open 
squares in Fig. [2]), we see that the separation of 
states in FLL is much larger than AMF (9 eV versus 
3 eV), with M = 7 the lowest energy for FLL but 
highest for AMF. This is a direct consequence of 
the magnetic penalty of AMF discussed previously. 
If / were increased above 1 eV (keeping the other 
parameters fixed), then AMF would begin to favor 
the M = 7 state by a small amount. 

Once J is turned on, the degeneracy is split, and 
the configurations with a particular M spread out 
around the J = value. The spread is especially 
large for the highly degenerate M = 1 value (from -5 
to 8 eV) , so that even if / were larger than the typical 
LSDA value (in which case, with J = AMF would 
favor a high spin state) the large spread of M = 1 
values would cause the low-spin states to be favored 
in AMF. This spread is entirely coming from the 
Ei term and is independent of the double-counting 
choice. Here we see for AMF a competition between 
J and /: J is actually preferring a low spin configu- 
ration, in contrast to the conventional wisdom that 
J increases the tendency for magnetism. We see this 
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same tendency occurs in FLL, as for J — the sepa- 
ration between M = 7 and M = 1 states is 9 eV, but 
with J = 1 this separation is reduced to 4 eV. Since 
in FLL the Hubbard U does not penalize magnetic 
states the way AMF does, the presence of J is not 
able to compete with /. This makes it clear why FLL 
is generally accepted to perform better for systems 
known to have high-spin states (e.g. Eu and Gd). 
Conversely, FLL may be less successful at modeling 
low-spin states. 

As mentioned previously, it is fairly common for 
theoretical studies to replace U and J with effective 
parameters U and J. For any double-counting term 
chosen, using these effective parameters will lower 
the energy of the high-spin state relative to the low 
spin state as compared to using U and J directly. 
With orbitals that are not highly localized, such as 
3d or 5/ states it may be the case with FLL that the 
reduction of the energy separation between high-spin 
and low-spin caused by using U and J would allow 
for significant competition between magnetism and 
kinetic energy in LSDA+U. 

We now have seen why and how FLL and AMF 
perform differently in assigning a magnetic moment. 
This may be of particular interest for studies of 
pressure-induced changes in magnetic moment, such 
as that seen in MnO^ without changes in orbital M 
occupancy. Applications of LSDA+U are more thor- 
oughly discussed in Sec. IVII Al 

Shown in Fig. [3] are scatter plots of the energies 
of all possible states for a given number of / elec- 
trons with integer occupations. SO is neglected, as 
it makes very minor changes to this picture by split- 
ting some degeneracies. The particle-hole symmetry 
of each functional is apparent. In Fl-nS and FLL, 
the ground state energy for N — 7 is roughly 3 eV 
lower than the next level, which are the (degenerate) 
ground states for N = 6 and 8. This is almost en- 
tirely due to the term depending on M 2 (either the 
J term in ([5]) or the Stoner term in FLL), because 
M is large and at its maximum with 7 spins aligned. 
In AMF low spin states can be seen at the low end 
of the range for configurations at each N; the high 
spin states for N—6 and 7 are disfavored by 6-7 eV. 
We see that the trend where AMF favors low-spin 
configurations and FLL favors high-spin configura- 
tions shown for N — 7 in Fig. [2] is present for all N. 
The large spread of values for low spin configurations 
(black and red circles) is seen clearly for AMF as 
they appear in both the lowest energy positions and 
the highest energy positions. The high-spin configu- 
rations (large open symbols and triangles) are in the 
middle of each distribution for N. For e.g. N = 5, 
counting from the lowest energy, M = 1 configu- 
rations are found first, followed by M — 3 config- 
urations then the M — 5 configurations are found 
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FIG. 3: (Color online) Scatter plot of all energies AE for all states in the (a) AMF(Fl-S), (b) FLL and (c) Fl-nS 
double-counting schemes, for U = 8, J = 1, and / = 0.75 (FLL and AMF only). Spin-orbit is neglected here. For 
AMF, low spin states (black and red circles) appear as lowest energy configurations for all N, but this is not the case 
for FLL or Fl-nS. The dashed lines indicate the mean energy over configuration for each N; note that the variation 
with N is much less for FLL than for the other two functionals. 



(with the trend reversing counting up to the highest 
energy states). In FLL, the lowest energy configu- 
rations for N + 1 7 are still the configurations with 
maximum spin for a given N, and states with lower 
spins are found in succession. Again using N = 5 as 
an example, the M = 5 configurations are lowest in 
energy, and then M = 3 configurations are seen at 
energies lower then M — 1 states. 



VII. DISCUSSION 

In this paper we have tried to clarify the behav- 
ior of the various functionals that are used in the 
LSDA+U method, we have compared the function- 
als formally in certain limits, we have presented the 
orbital potentials that arise, and we have analyzed 
the total energy corrections that LSDA+U function- 
als apply to LSDA total energies, given a set of oc- 
cupation numbers. The Fl-nS functional which was 
originally introduced strongly favors spin-polarized 
states as does the commonly used FLL functional. 
The other most commonly used functional besides 
FLL, Fl-S (AMF), has characteristics that tend to 
suppress moment formation or reduce the magnitude 
of the moment. When analyzed, this AMF func- 
tional shows positive energy penalties to magnetism 
that compete with the magnetic tendencies of the 
LSDA functional, and when J > non-magnetic 
solutions become even more likely to win out. We 
have provided a short analysis of the behavior when 
J = is used. While this case is instructive, we ad- 
vise against its use; it is just as simple to do the full 
J ^ calculation. 



When LSDA+U is applied to correlated insula- 
tors in the strong coupling regime, it provides a very 
good picture of the system at the band structure (ef- 
fective one-electron) level. The initial successes in- 
clude the 3d transition metal monoxides MnO, FeO, 
CoO, and NiO, for which the LSDA description is 
very poor. Other early successes included the insu- 
lating phases of the layered cuprates that become 
high temperature superconductors when doped, and 
the unusual magnetic insulator KCUF3, which was 
the first case where crucial orbital ordering was re- 
produced. LSDA+U is not a satisfactory theory of 
single particle excitations of such systems, but never- 
theless provides a realistic picture of the underlying 
electronic structure. 



The more interesting, and more difficult, cases 
now lie between the strongly correlated limit of 
wide-gap magnetic insulators and weakly correlated 
regime that is well described by LSDA. Some of these 
are metals, some are unconventional insulators, and 
many lie near the metal-insulator borderline. It is 
for these intermediate cases that it becomes essen- 
tial, if applying the LSDA+U approach, to under- 
stand what the method is likely to do, and especially 
to understand the tendencies of the various choices 
of functional. This is what we have tried to clar- 
ify in this paper. As a summary, we will provide 
an overview of an assortment of results that have 
appeared in the literature for systems that lie some- 
where in the intermediate correlation regime. 



A. Examples of LSDA+U behavior from 
applications 

1. Strongly correlated insulators. 

Cuprates. The insulating phase of the cuprate 
class of high temperature superconductors com- 
prised the "killer app" that served to populariz o 1 ' 20 
the LSDA+U method, and in the intervening years 
the method has been applied to cuprates and other 
correlated insulators too many times to cite. Sim- 
ply put, in cuprates it produces the Cu d 9 ion and 
accompanying insulating band structure. 2 ^ 21 The 
hole resides in the d x 2_ y 2 orbital and is strongly hy- 
bridized with the planar oxygen p a orbitals, as much 
experimental data was indicating. 

MnO. Experimentally, MnO shows at room tem- 
perature a moment collapse from M = 5 to M = 1 
(or less), a volume collapse, and an insulator- to- 
metal transition, near 100 GPa; this is the classic 
Mott transition. Within LSDA, the moment de- 
creases continuously with decreasing volume^ 2 - from 
the high spin (HS) state to a low spin (LS) state. 
The insulator-to-metal transition occurs at much too 
low a pressure (without any other change). A vol- 
ume collapse is predicted, although the pressure is 
significantly overestimated (150 GPa). 

The application of LSDA+U in its FLL flavor 
has been applied and analyzed in detail^ and pro- 
vides a different picture in several ways. The am- 
bient pressure band gap is improved compared to 
experiment. The volume collapse transition occurs 
around 120 GPa and is accompanied by a moment 
collapse from M = 5 to M = 1. The nature of 
this (zero temperature) transition is insulator to 
insulator, while the experimental data indicate an 
insulator-to-metal transition at room temperature. 
The zero temperature transition might indeed be 
insulator-to-insulator; such a phase transition would 
be a type that LSDA+U should work well for. It 
is also possible that the static mean-field approxi- 
mation underlying LSDA+U, which favors integer 
occupations and hence insulating solutions, has too 
strong a tendency and fails to describe this tran- 
sition. This question could be settled by studying 
experimentally the Mott transition at low tempera- 
ture 

Even more unexpected than the insulator to insu- 
lator aspect is the LSDA+U prediction is that the 
low spin state has an unanticipated orbital occupa- 
tion pattern^ being one in which every 3d orbital 
remains singly occupied (as in the high spin state), 
but spin in two orbitals antialign with those in the 
other three orbitals. This state is obtained simply 
from the M = 5 high spin state by flipping the 
spins of two of the orbitals. The resulting density 
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remains spherical, but the spin density exhibits an 
angular nodal structure leading at the same time to 
a high degree of polarization of the spin-density but 
a low total moment (M = 1). This solution (be- 
ing the high pressure ground state in LSDA+U) can 
be traced^ back to the interplay between symmetry 
lowering due to the antiferromagnetic order (cubic 
lowered to rhombohcdral) and the anisotropy part 
of the interaction Eq. 1201 The symmetry lowering 
lifts the cubic grouping (t-2 g and e g manifolds), thus 
allowing a higher number of allowed occupation pat- 
terns. 

The anisotropic part of the interaction is 
responsible^ for Hund's second rule ordering of 
states, which has the tendency to increase the mu- 
tual distance of each pair of electrons. If the over- 
all energetics (band broadening and kinetic effects) 
reduce the gain of energy due to spin-polarization, 
then Hund's first rule may become suppressed and 
the result is a low spin state. The anisotropic inter- 
action is however not influenced by this suppression, 
since it is a local term proportional to a parameter 
J. It will enforce a Hund's second rule like separa- 
tion of the electrons under the low spin condition, 
and thus can be shown to result exactly in the oc- 
cupation pattern observed for MnO. In a sense the 
low spin state is an example of Hund's second rule 
without Hund's first rule. 

FeO, CoO, MO. Together with MnO, these 
classic Mott (or 'charge transfer') insulators 
have been prime applications of the LSDA+U 
methodi^^ii^i 2 ^ The behavior of the open 3d shell 
in these compounds has not been analyzed in the 
detail that was done for MnO, however. 



2. Metals 

Correlated metals involve carriers that can move, 
hence they invariably involve fluctuations, in occu- 
pation number, in magnetic moment, in orbital oc- 
cupation, etc. It cannot be expected that a self- 
consistent mean field treatment such as LSDA+U 
can answer many of the questions raised by their 
behavior. However, there is still the question of 
whether LSDA+U can provide a more reasonable 
starting point than LSDA alone in understanding 
these metals. In our opinion, this remains an open 
question, but one for which some evidence is avail- 
able. 

The Fe-Al system has provided one platform for 
the application of LSDA+U to moderately corre- 
lated metals. The systems treated include the Fe 
impurity in Al (Kondo system, experimentally) , and 
the compounds F83AI, FeAl, and FeAl 3 . The calcu- 
lated behavior is too complex to summarize here. 



The LSDA+U result will, generally speaking, be 
likely to give a good picture of a Kondo ion when 
it produces an integer-valent ion with a large value 
of U. Both FLL and AMF functionals have been ap- 
plied in this regime ) 27 ' 28 with substantially differ- 
ing results, leading one to question whether either 
is more realistic than simple LSDA. Results are also 
sensitive to volume, i.e. whether using the exper- 
imental lattice constant or the calculated equilib- 
rium value, and the calculated equilibrium is differ- 
ent from LSDA and LSDA+U. One result was that, 
for moderate LV e ~ 3-4 eV, AMF strongly reduces 
the magnetic moment, while FLL does notj^ 8 - An- 
other application found that the magnetism disap- 
peared within a certain range of intermediate values 
of Up e , that is, it was magnetic around small Up e 
and also again at large coupling, 27 but non-magnetic 
between. 



3. Moderately strongly interacting oxides. 

Trying to address seriously the electronic struc- 
ture of intermediate coupling oxides, which are of- 
ten near the metal-insulator transition, is a challenge 
that has begun to be addressed more directly. The 
peculiar Na^OoC^ system, which becomes supercon- 
ducting when hydrated (water intercalates between 
C0O2 layers) is one example. One set of studies 
showed no appreciable difference between FLL and 
AMF, 29 with both predicting charge disproportion- 
ation on the Co ion for a?=| and | for U ~ 2.5-3 
eV. It is likely that this compound presents a case 
where the interplay between LSDA and U has ef- 
fects that are not fully understood. Also, it is un- 
clear why there is so little difference between the 
FLL and AMF functionals in this system. 

The compound Sr2Co04 is another example. 
Both functionals show a collapse of the moment 30 
around U — 2.5 eV, related to the metal-half metal 
transition that occurs, but the result for the mo- 
ments (M(AMF) < M(FLL)) bears out the ten- 
dency of AMF to penalize magnetic moments.. The 
fixed spin-moment calculation in Fig. 9 in Ref. [10 
is instructive too, showing the competition between 
LSDA magnetic energy and AMF magnetic penalty. 
Also it shows the creation of local minima around M 
= integer values that LDA+U introduces. 



4- f electron materials. 

4f systems. These metals often display the cor- 
related electron physics of a magnetic insulator at 
the band structure level: background conduction 
bands provide the metallic nature, while the corre- 
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lated states have integer occupation. The LSDA+U 
method seems to be a realistic method for plac- 
ing the / states closer to where they belong (away 
from the Fermi level). Gd is a good example, which 
has been studied at ambient pressure and compared 
to photoemission data 24 and magnetic dichroism 
datai^L 3 ^ The LSDA+U method has also been ap- 
plied up to extremely high pressure to assess where 
the 'Mott transition' in the 4/ bands is likely to oc- 
cur. The LSDA+U method has also been applied 
to heavy fermion metals, for example Cu and U 
compounds^ 3 - PrOs 2 Sbi2, [13 andYbRh 2 Si 2 H|. In 
such systems the LSDA+U method may even pro- 
vide a good estimate of which itinerant states at 
the Fermi level are strongly coupled to the localized 
/ states, i.e. the Kondo coupling matrix elements. 
These 4/ systems may become heavy fermion metals 
(YbRh 2 Si 2 ) or novel heavy fermion superconductors 
(YbAlB^, or they may remain magnetic but other- 
wise rather uninteresting metals (Gd). 

5/ systems. A variety of application of the 
LSDA+U method to 5/ systems, and especially Pu, 
have been presented ! 36 ' 37 ' 38 ' 39 ' 40 Given the complex- 
ity of the phase diagram of elemental Pu, together 
with claims that dynamic correlation effects must be 
included for any realistic description of Pu, a more 
critical study of Pu would be useful. 
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APPENDIX A: CALCULATION OF THE 
STONER I FOR 3d AND 4/ SHELLS 

The Stoner parameter / is a well established quan- 
tity. For metals its value is obtained by a second 
order expansion of the LSDA xc-energy around the 
non-magnetic solution, resulting in a Fermi surface 
averaged integral of the radial wave functions with 
the xc-kernel4i LSDA+U is usually applied to de- 
scribe insulating states, where the Fermi surface van- 
ishes. In the context of discussing the LSDA contri- 
bution to the energy of a correlated d- or /-shell it 
is more natural to consider the energy contribution 
from the localized shell. This leads to a derivation 
of the Stoner-/ similar to the formulation of Janak 
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FIG. 4: Shell-Stoner integrals for the 3d and 4/ atoms. 
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but adapted to atom- like situations. 

Seo presented^ the second order perturbation 
theory of the spin polarization in DFT, which results 
in explicit expressions for the shell exchange param- 
eter I n i that are applicable to atom like situations. 
In this work a numerical estimate for I n i was derived 
indirectly from exchange splittings and spin polar- 
ization energies taken from DFT calculations. The 
idea behind this perturbation theory, the expansion 
the xc-energy around the spherically averaged non- 
magnetic density of the shell under consideration, 
was also discussed in the appendix of Kasinathan et 



al^ and leads to AE XC 
Stoner integral 
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The last expression applies for a local approximation 
(viz. LSDA) to E xc . Ko^r 1 ) is a magnetization- 
magnetization interaction, directly analogous to the 
second functional derivative of the DFT potential 
energy with respect to n(r), which is the Coulomb 



interaction e /|f 
from E xr . 



plus an 'xc interaction' arising 



(A2) 



For a more detailed discussion of the parameter 
Ini we performed LSDA calculations for free atoms 
and ions and explicitly calculated from Eq. (|A1[) . 
It turns out that AE XC (M) given above is by far the 
largest M-dependent term of the energy expansion. 
The spin polarization energy of isolated atoms/ions 
with spherical M is well described by this estimate 
with an error smaller than 5 — 10%. The result- 
ing shell-Stoner integrals I n i have very similar values 
compared to the ones obtained from the theory for 
the metallic situation. (Note, however, that there is 
a factor of 2 difference in the definition of the Stoner 
I in some of the publications.) 



For the 3d transition element series we get values 
I n i ranging from 0.62 eV for Sc to 0.95 eV for Zn (see 
Fig. 4). These values increase across the series by 
« 0.15 - 0.20 eV, when the exchange only LSDA is 
used, pointing to a reduction due to (LDA-type) cor- 
relation effects when the full xc-kernel is used. For 
the 4/ series the shell-Stoner integrals vary from 0.58 
cV for Ce to 0.75 eV for Yb. The LDA correlation 
effects amount to 10% of these values. The values 
obtained depend on the choice of the reference sys- 
tem, which serves as zeroth order in the functional 
expansion. For instance for the 3 + -ions of the 4/- 
series J4/ is increased by 6 — 20% with respect to the 
neutral atoms. 
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